Lithium Adsorption on Zigzag Graphene Nanoribbons 



Chananate Uthaisar, Veronica BaroneO and Juan E. Peralta 
Department of Physics, Central Michigan University 
Mount Pleasant, MI 48859. 
(Dated: October 27, 2009) 

We have studied the adsorption of Li atoms at the hollow sites of graphene nanoribbons (zigzag 
and armchair), graphene, and fullerenes by means of density functional theory calculations including 
local and semilocal functionals. The binding energy of a Li atom on armchair nanoribbons (of about 
1.70 eV for LSDA and 1.20 eV for PBE) is comparable to the corresponding value in graphene (1.55 
and 1.04 eV for LSDA and PBE, respectively). Notably, the interaction between Li and zigzag 
nanoribbons is much stronger. The binding energy of Li at the edges of zigzag nanoribbons is about 
50% stronger than in graphene for the functionals studied here. While the charge transfer between 
the Li adatom and the zigzag nanoribbon significantly affects the magnetic properties of the latter 
providing an additional interaction mechanism that is not present in two-dimensional graphene 
or armchair nanoribbons, we find that the morphology of the edges, rather than magnetism, is 
responsible for the enhanced Li-nanoribbon interaction. 

PACS numbers: Valid PACS appear here 



I. INTRODUCTION 



The recent experimental realization of atomically thin, 
long, and narrow strips of graphene (graphene nanorib- 
bons, GNRs) [E H H has sparked an intense research 
effort toward the understanding of these novel materi- 
als. Experimental evidence of ballistic electronic trans- 
port, large phase coherence lengths, and current den- 
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ior, [13lll4 Ul5lll6l ll7ll and electronic structure engineer- 
ing capabilities [H M, SH S3, HI H H identify 
low-dimensional graphene as one of the most promising 
materials for novel electronic and spintronic devices. Un- 
like two-dimensional graphene, GNRs present electronic 
confinement in the transverse direction giving rise to pe- 
culiar electronic properties, similar to the case of sin- 
gle walled carbon nanotubes (SWNTs). As in SWNTs, 
these properties depend strongly on the crystallographic 
orientation of their main axis, i.e. armchair or zigzag. 
In addition to their remarkable electronic and magnetic 
properties, GNRs present a distinctive chemistry due to 
the presence of reactive edges. [H H|| [2tJ 

The interaction between Li atoms and carbon-based 
materials has been a subject of intensive study moti- 
vated by their potential for more efficient Li-ion batteries 
and hydrogen storage media. While the detailed interac- 
tion mechanism is still controversial, [281 [29( it is com- 
monly accepted that it presents a mostly ionic character 
with a substantial electronic charge transfer from the Li 
atom to the graphitic surface. [3(J, HH, HH, [33| In a re- 
cent work, Yoo et al. presented experimental evidence of 
a higher Li storage capacity in graphene with respect to 
graphite. [3(| Li adsorption on fullerenes has been studied 
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theoretically by Sun et al. [32j In that work the authors 
found that Lii 2 C6o can store a high density of molecular 
hydrogen. Recently, Ataca et al. presented a theoretical 
study of the adsorption of Li atoms on graphene with no- 
tably large Li densities. [34| However, the binding energies 
of Li atoms of such systems, with high Li coverage, are 
significantly lower than that of a less dense Li-graphcnc 
system due to the electrostatic repulsion between Li ions. 

In view of all the recent focus on Li storage in low- 
dimensional carbon materials, we present this work with 
the aim of studying the adsorption mechanisms and 
binding energies of Li atoms on graphene nanoribbons 
and assess their Li intake capacity compared to other 
low-dimensional carbon materials. In what follows, we 
present theoretical evidence, based on density functional 
theory calculations, of a stronger interaction between Li 
atoms and zigzag nanoribbons in comparison with cither 
graphene or armchair nanoribbons. 



II. METHODOLOGY 

We have carried out all calculations utilizing a devel- 
opment version of the Gaussian suite of programs. 35] In 
this program, solid state calculations are performed using 
all-electron Gaussian basis sets and periodic boundary 
conditions in one, two, or three dimensions. This flexi- 
bility allows for the treatment of low-dimensional struc- 
tures, such as graphene or graphene nanoribbons, avoid- 
ing replicas in the z direction. Calculations have been 
performed using an spin-polarized Kohn-Sham approach 
within the local spin density approximation (LSDA) 
[3(| [3^ | and the generalized gradient approximation of 
Perdew et al. [MlHl We chose the double-zetha 6-31G** 
basis set. [4(| Despite the fact that the ionic character 
of the Li-grpahene interaction should be qualitatively 
well reproduced by either LSDA or GGA functionals, 
we studied the case of Li adsorption on zigzag nanorib- 
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FIG. 1: Schematic diagrams of the supercells used for a) two- 
dimensional graphene, b) armchair graphene nanoribbon, and 
c) zigzag graphene nanoribbon. 



bons using an hybrid density functional. As upon Li 
adsorption the systems turn into metals, conventional 
hybrid functionals present convergence problems due to 
the portion of Hartree-Fock type of exchange included in 
the exchange-correlation energy. An efficient alternative 
to deal with this problem is the screened-exchange hy- 
brid functional, HSE, developed by Heyd, Scuseria, and 
Ernzerhof . El, 

The reciprocal space integration has been performed 
in a 10x10 uniform k-point grid for graphene and in a 
35 uniform k-point grid for the one-dimensional nanorib- 
bons. All the structures have been fully relaxed to ac- 
count for any rearrangement in the graphitic surfaces 
due to the presence of the Li atom until the maxi- 
mum and root mean square atomic forces are less than 
0.02 eV/Aand 0.015 eV/A, respectively, and the max- 
imum and root mean square atomic displacements be- 
tween consecutive iterations are less than 10~ 3 Aand 
6 10~ 4 A, respectively. 

We have studied Li adsorption on three representa- 
tive graphene surfaces (shown in Figure 1). The unit 
cell for the first system consists of an isolated (4x4) su- 
percell of two-dimensional graphene (without replica in 
the z direction) with 32 C atoms and one Li atom with 
a Li. ..Li distance between cells of 0.98 nm (Figure 1-a). 



The interaction Li-graphene in this particular configu- 
ration has been studied before and therefore provides a 
good reference system for comparison purposes. 0, EH] 
The second system consists of an isolated supercell of 
an armchair nanoribbon 1.58 nm wide, containing 56 C 
atoms, 8 H atoms passivating the edges, and one Li atom 
(Figure 1-b). Within this configuration the Li.... Li dis- 
tance between cells is about 0.86 nm. Finally we have 
chosen an isolated supercell of a zigzag nanoribbon 1.57 
nm wide, containing 64 C atoms, 8 H atoms passivating 
the edges and one Li atom in such a way that the Li.... Li 
distance between cells is 0.99 nm (Figure 1-c). For com- 
parison purposes, we have also considered the case of 
Li adsorption on fullcrcncs, which presents one of the 
strongest adsorptions reported between Li and graphitic 
materials. [32j 

It has been shown in several studies that alkali metal 
adsorption on graphitic surfaces takes place preferentially 
on top of the hexagon (hollow site) instead of on top of 
a C atom or on top of a C-C bond (bridge site). 0, 
Therefore, we have studied the interaction between a Li 
atom on the hollow site of the graphitic materials in the 
three configurations described above and checked that 
the hollow sites are more stable for the case of the zigzag 
nanoribbon (Figure 1-c). 

Since we are employing a finite localized basis set, our 
calculations are subject to basis set superposition error 
(BSSE). Although the correction for BSSE can be de- 
manding due to the large number of configurations stud- 
ied here and the different spin polarized ground states of 
the constituents, we can estimate its magnitude by calcu- 
lating it using the counterpoise method [4 7j in the case of 
Li adsorbed on two-dimensional graphene. In this case, 
basis set truncation introduces an error in the computed 
binding energy smaller than 4%. For this reason, we ne- 
glect the BSSE correction in the rest of this work. 



III. RESULTS 

1. Graphene, fullerenes, and armchair nanoribbons 

Several papers appeared in the literature present- 
ing DFT calculations in (4x4) graphene containing one 
adatom per 32 C atoms. For instance, Khantha at al.0] 
report a binding energy (Eb) of 1.60 eV; Valencia et 
al. 0| report the Li binding energy in this system to 
be 1.55 eV with LSDA and 1.01 eV with the general- 
ized gradient approximation functional of Perdew-Burkc- 
Ernzerhof, PBE, while Chan et al. report a binding en- 
ergy of 1.10 eV with PBE as well. [46j More recently, 
Ataca et al. calculated the binding ener gy for Li on 
(4x4) graphene to be 1.93 eV with LSDA.Q For this 
system, we obtain a binding energy of 1.55 eV (LSDA) 
and 1.04 eV (PBE) in good agreement with the values 
reported by Valencia et al. 45 1, Khantha at al.,0] and 
Chan et al.0 

The interaction of Li with Cqq is stronger than with 
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TABLE I: Calculated Li-hexagon distances, magnetic moments per cell and Mulliken atomic charges on the Li atom adsorbed 
on graphene, Ceo, and the four possible configurations of Li on the hollow site of a zigzag GNRs. 
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FIG. 2: LSDA and PBE Li binding energies on the hollow 
sites of the different systems considered in this work: (4x4) 
graphene, Ceo, armchair nanoribbons with 1 < i < 6, and 
zigzag nanoribbons with 1 < i < 4. 



armchair GNRs. Both hollow sites (on top of the pen- 
tagon and on top of the hexagon) are degenerate with 
a binding energy of 1.94 eV for LSDA and 1.42 eV for 
PBE. This value, calculated using PBE, is smaller than 
the 1.80 eV obtained by Sun et al. with the PBE func- 
tional. 

Due to the quasi-one-dimensional nature of the 
nanoribbons, there are several possible distinctive con- 
figurations for Li on the hollow sites. These different 
possibilities are schematized in Figures 1-b for armchair 
and 1-c for zigzag nanoribbons and denoted hex — i with 
1 < i < 6 for armchairs and 1 < i < 4 for zigzags. 
As shown in Figure 2, the interaction of Li with arm- 
chair nanoribbons is to some extent stronger than with 
graphene but still weaker than with Cgo ■ Even when the 
first neighbor Li... Li distance in the armchair nanoribbon 
chosen here is slightly smaller than in the 4x4 graphene 
supercell, the Coulomb repulsion between ions is stronger 
in graphene. Li atoms in graphene have four first neigh- 
bors (two-dimensional structure) while nanoribbons only 
have one first-neighbor (one-dimensional structure) . This 



difference increases the electrostatic repulsion between 
ions in graphene. As in the case of graphene, Li ad- 
sorption on armchair nanoribbons does not induce spin 
polarization in these systems. 

There is a slightly stronger interaction between Li and 
armchair GNRs when the Li atom is on the hex-1 position 
that corresponds to the hexagons at the edges (Figure 2). 
In all other cases, there are no significant differences in 
the binding energies for the width studied here. For wider 
ribbons, we expect Ef, to tend asymptotically to a value 
somewhat larger than the one calculated in graphene due 
to their one-dimensional nature and its consequence on 
the number of first-neighbor Li atoms. 



2. Zigzag nanoribbons 

One of the most remarkable results found here is that 
the interaction of Li with zigzag GNRs is much stronger 
than with armchairs. As shown in Figure 2, the binding 
energy of Li at the edges of zigzag nanoribbons is more 
than 50% stronger than in graphene for both functionals, 
LSDA and PBE. This interaction slowly weakens as the 
adsorption position is shifted toward the center of the 
nanoribbon suggesting that narrow zigzag GNRs will be 
the most ideal to attain larger ion concentrations. While 
graphene and armchair nanoribbons do not develop mag- 
netism upon Li adsorption, the charge transfer between 
the Li atom and the zigzag nanoribbon affects signifi- 
cantly its magnetic properties providing yet another in- 
teraction mechanism that is not present in graphene and 
armchair nanoribbons. 

It is worth mentioning the effect of geometry relax- 
ation in the binding energy for the different adsorption 
sites. We find that this effect is rather small across the 
hex — 1, hex — 2, hex — 3, and hex — 4 series and accounts 
for, at most, 0.06 eV of the total binding energy, calcu- 
lated as the difference between the binding energies of 
the fully relaxed system and the rigid model approach at 
the equilibrium distance. This indicates that the stronger 
binding energy at the edges does not originate in defor- 
mations of the carbon backbone. 

To further investigate the interaction of Li and zigzag 
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nanoribbons we calculated the binding energies for the 
hex — i series with the hybrid HSE functional. We obtain 
binding energies of 1.55, 1.27, 1.16, and 1.13 eV for hex—i 
with z=l,2,3,4, respectively. These results confirm the 
trend predicted by the LSDA and PBE functionals. 

Zigzag edges exemplify the importance of edge effects 
in honeycomb lattices. [26l . |49| Zigzag nanoribbons are 
expected to present a spin-polarized ground state charac- 
terized by an antiferromagnetic spin arrangement (AFM) 
with opposite spins at each edge. (5(| The high spin state 
solution, with all spins ferromagnetically aligned (FM), 
is higher in energy than the AFM state by 10 meV/edge 
atom for a 1.8 nm wide ribbon. This magnetic behav- 
ior has been predicted to present robustness with respect 
to some edge defects and impurities. [l(| [III However, 
we find that when Li is adsorbed on the zigzag nanorib- 
bon, the charge that is transferred from the Li atom to 
the nanoribbon quenches the magnetization of the car- 
bon atoms in the vicinity of the adsorption site in such a 
way that magnetization at both edges is no longer com- 
pensated, resulting in a net magnetic moment per cell 
that strongly depends on the adsorption position. Inter- 
estingly, except in the hex — 1 case that quenches the 
magnetization at the edge in which the Li atom is ad- 
sorbed, the AFM and FM solutions become almost de- 
generated [5l[ while presenting significantly different to- 
tal magnetic moments as shown in Table 1. These effects 
are shown in Figure 3 where the LSDA spin density maps 
0.05 A above the surface of the AFM solution (ground 
state) of the pristine graphene nanoribbon and the dif- 
ferent positions of the Li atom for the Li-GNR systems 
are presented for every adsorption site. We show the FM 
solution for the Li-GNR system as the AFM solution is 
very similar to the pristine ribbon but with a very small 
spin density in the vicinity of the adsorption site, just as 
in the FM case. In all cases, the Li adatom suppresses 
the spin polarization surrounding it, in such a way that 
spin-polarized edges do not appear for the hex-1 position 
on the Li side but magnetism slowly develops for hex-2 
and increases for hex-3 and hex~4- It is worth mention- 
ing at this point that a spin compensated calculation us- 
ing the local density approximation produces a nonmag- 
netic metallic solution. This solution is higher in energy 
than the spin-polarized solution by about 70 meV/cell 
for hex — 1 and 108 meV/cell for hex — 4. In the pris- 
tine case, the energy difference between the nonmagnetic 
and AFM solution is about 160 meV/cell. Therefore, if 
we do not consider spin-polarization at all and perform 
all the calculations but the isolated Li atom using the 
spin-restricted LDA functional, the binding energy of the 
Li adatom in the different positions remains almost the 
same as in the spin-polarized case. These results indi- 
cate that edge morphology rather than magnetization is 
responsible for the enhanced interaction between Li and 
zigzag nanoribbons and a that a stronger binding is ex- 
pected in zigzag nanoribbons regardless their magnetic 
nature. 

In order to understand the difference in the adsorption 
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FIG. 3: LSDA spin density maps 0.05 A above the surface 
of the zigzag graphene nanoribbon for different positions of 
the Li atom. Blue and red colors represent opposite spin 
polarization. 



strength of Li in zigzag and armchair nanoribbons, we 
analyze the total and partial density of states (DOS) in 
the isolated systems. In Figure 4 (b) we present the total 
DOS of a zigzag and armchair nanoribbons. The DOS 
corresponding to the zigzag nanoribbons is much larger 
for states close to the Fermi level (Ep) than in armchair 
ribbons. This larger number of states makes the charge 
transfer from the Li atom more accessible in the case of 
zigzag than in armchairs GNRs. Moreover, the stronger 
interaction at the edges of zigzag nanoribbons can also 
be understood by considering the partial DOS. In Figure 
4 (c) and (d) we show the total and partial DOS for arm- 
chair (c) and zigzag (d) nanoribbons. The partial DOS is 
constructed by only considering p orbitals perpendicular 
to the GNR from C atoms belonging to the different lay- 
ers shown in Figure 4 (a). In this way, we can separate 
the contributions to the total DOS from the edges (layer 
1) and the center part of the ribbon (layer 4). As seen 
in Figure 4 (c) and (d), while for the armchair GNR all 
layers contribute about the same to the total DOS in the 
vicinity of the Fermi level, the corresponding contribu- 
tions in zigzag nanoribbons are significantly different. In 
the zigzag case, both, valence and conduction bands close 
to Ep are dominated by edge states (from layer 1) and 
layer 2 while layers 3 and 4 present significant smaller 
contributions. This explains why the binding energy of 
the adatom is much stronger at the edges and slowly de- 
creases toward the center of the zigzag nanoribbon. 



IV. CONCLUSIONS 

In summary, we have studied the adsorption of Li 
atoms at the hollow sites of fullerenes, graphene, and 
graphene nanoribbons (zigzag and armchair) by means 
of density functional theory within the local spin den- 
sity and generalized gradient approximations. Li in- 
teracts with armchair nanoribbons and two-dimensional 
graphene through the same charge transfer mechanism, 
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with binding energies per adatom of about 1.70 eV and 
1.55 eV (LSDA) and 1.04 eV and 1.20 eV (PBE), respec- 
tively. Li interacts with zigzag nanoribbons in a much 
stronger way. The binding energy of Li at the edges 
of zigzag nanoribbons is 2.27 eV (LSDA) and 1.70 eV 
(PBE), more than 50% stronger than in graphene. The 
binding energy progressively decreases as the adsorption 
position is shifted toward the center of the nanoribbon 
suggesting that narrow zigzag GNRs will be most ideal 
to attain larger ion concentrations. 

While the charge transfer between the Li adatom and 
the zigzag nanoribbon affects significantly the magnetic 
properties of the latter providing an additional interac- 
tion mechanism that is not present in graphene or arm- 
chair nanoribbons, we find that the morphology of the 
edges, rather than magnetization, is responsible for the 
enhanced Li adsorption. These results illustrate the im- 
portance of controlling the edges of GNRs with atomic 
precision in order to maximize their potential for tech- 
nological applications. This precision has been demon- 
strated to be experimentally achievable.^ 
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FIG. 4: (a) Schematic representation of the different C atom 
layers in zigzag and armchair nanoribbons. (b) Total DOS for 
the armchair and zigzag nanoribbon studied in this work, (c) 
Total and partial DOS for the armchair GNR with contribu- 
tions from the 4 different layers shown in (a), (d) Total and 
partial DOS for the zigzag GNR with contributions from the 
4 different layers shown in (a). 
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